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Abstract. I study hard collisions between unpolarized protons and antiprotons where a lepton-antilepton 
pair is detected in coincidence with a final proton-antiproton pair, and no more particles are produced, 
in the regime 10 GeV 2 < s < 1000 GeV 2 , M > 4 GeV, q T < 3 GeV/c. The present work is centered 
on azimuthal asymmetries. Because of momentum conservation, a Boer-Mulders term in the momentum 
distribution of a quark implies a balancing effect in the momentum distribution of some spectators. This 
produces azimuthal asymmetries of the final hadrons. To analyze this, I have organized a parton-level 
MonteCarlo generator where a standard cos(20)-asymmetry of the dilepton distribution is produced, thanks 
to a soft rescattering process between an active quark coming from a hadron and a spectator anti-diquark 
coming from the other hadron. This produces cos(2</>)-asymmetries of the final hadron pair. Hadron and 
lepton asymmetries have the same size. 

PACS. 13.85.Qk Drell-Yan - 13.85.Hd Inelastic processes with multi-particle final states - 13.88.+e Po- 
larization in interaction and scattering 



1 Introduction 

The rising interest for the so-called spin-physics in inclu- 
, sive hadronic reactions (see e.g. [T] or [5]) has among its 
consequences the planning of Drell-Yan experiments [3J|3] 
at intermediate energies (squared cm. energy s <C 1000 
GeV 2 ). One is the PANDA experiment [5] where protons 
and antiprotons collide with variable s < 30 GeV 2 (see [5] 
, for an extensive description of the experiment goals and 
methods, and [7] for a shorter summary). This experiment 
has a broad spectrum of objectives. These include unpo- 
larized Drell-Yan, finalized at a precision measurement of 
' the dilepton cos ( 2 </>)— asymmetry for s w 30 GeV 2 , and 
dilepton mass M > 2 GeV/c 2 . 

The cos (20)— asymmetry, or v— asymmetry, was first 
discussed in [8j|9]. It was found nonzero, with size 5-20 %, 
in 7r~— nucleus experiments [TUIITTllT2lll3l[Tl] . where s was 
between 100 and 600 GeV 2 and both x and x ~ 0.1. It 
was recently found small or zero in the pp experiment [15] 
at much larger s and much smaller average x, x. PANDA 
will measure it in a peculiar large— x regime, with large 
statistics. Both muon pairs and electron pairs will be de- 
tected. 

Apart for specific kinematics, another peculiarity dis- 
tinguishes PANDA from all the previous fixed-target Drell- 
Yan experiments. Because of the multi-purpose nature and 
design of this apparatus, it will detect the class and mo- 
mentum of almost all of the reaction products. 



In the past, the most appealing feature of the Drell- 
Yan measurement (in the case of dimuon detection) has 
just been the need not to care other particles apart for the 
pair. With a thick layer of hadron-absorbing mate- 
rial and muon spectrometers downstream, the structure 
of a fixed-target Drell-Yan apparatus is relatively simple 
(see [I61IT7] for reviews). For this reason, we have very 
little information about what is produced in Drell-Yan 
measurements, in coincidence with the dilepton pair. The 
author is only aware of charge multiplicity measurements 
at ISRfTMTg]. at the large y/s ~ 100 GeV (where a large 
number of particles was produced, of course). 

To analyze what PANDA could find concerning Drell- 
Yan fragments, in a previous work 20 I presented a sim- 
ulation of Drell-Yan events at s = 30 GeV 2 , with min- 
imum dilepton mass 2 GeV/c 2 and minimum transverse 
momentum 0.8 GeV/c. That simulation was performed 
with Pythia-8[21 . I found some surprising results: 

1) Almost all the events contain one (and only one) 
nucleon-antinucleon pair. 

2) 50 % of the events only contains a nucleon-antinucleon 
pair (apart for the dilepton, of course). These pairs are 
equally divided into pp and nfi pairs. 

3) In most of the other events the NN pair is accom- 
panied by not more than two light particles (charged pions 
or photons, possibly from ir° decay). 

4) An insertion of the simulated events into the PANDA 
acceptance shows that in the case of "dilepton + dipro- 
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ton" events with a detected dilepton, half of the pp pairs 
will be fully detected. This is about 10 % of the total rate 
of detected Drell-Yan events in PANDA conditions. 

5) In another 10 % of detected Drell-Yan events we will 
only detect a single proton (often) or a single antiproton 
(more rarely), but we will be able to identify a "dilepton 
+ diproton" exclusive event. 

This opens the possibility to analyze 20 % of the Drell- 
Yan events as exclusive events, with some interesting as- 
sociated observables and with the perspective of accessible 
modeling (because of the small number of involved parti- 
cles). This work will be centered on an analysis of the 
"dilepton + diproton" events, including azimuthal asym- 
metry effects. Clearly, the standard properties of these 
events may be simulated via Pythia or some other known 
multi-purpose code as I have already done in [5D]. Here 
my primary goal is to analyze those unknown properties 
of the proton-antiproton pair, associated with the dilep- 
ton azimuthal cos (20)— asymmetry. To this aim, I have 
organized a specific parton-level generator code. 

The extra detected proton-antiproton pair enriches the 
analysis with two more vectors: pp and pp. These may be 
combined into the two pairs pp ±pp. In 4-particle events 
only the relative momentum pp — pp is independent of 
the dilepton momentum q (that is equal to the sum of 
the lepton momenta and to minus the sum of the frag- 
ment momenta). So, here I will focus on the observables 
related with the correlations among the fragment relative 
momentum, the lepton relative momentum and q. 

Although my main interest is towards PANDA, I will 
consider a hypothetical experiment with s = 100 GeV 2 
instead of 30 GeV 2 as in PANDA (a possibility in this 
kinematic range could be Drell-Yan at COMPASS [2"2"]). 
This will allow me to use several simplifying approxima- 
tions, based on the fact that the longitudinal momenta of 
the hadron constituents are for most events much larger 
than the transverse momenta. In the PANDA case this 
is not true, and the analysis would require nonstandard 
devices, and infrared parameters. 

1.1 Notations 

I will use the letter "p" to indicate observable momenta 
(leptons, proton, antiproton) and the letter "fc" for the 
momenta of the partons (quark, diquark, antiquark, an- 
tidiquark). In detail, p M and p^ indicate lepton momenta, 
kft and fc M (anti)quark momenta, while for the other par- 
ticles a suffix is present. 

The vector q M is the sum of the lepton momenta, and 
also of the quark momenta, and minus the sum of the 
spectator momenta and of the final hadron momenta. I 
use both q = y'q^q^ and M to indicate the same thing, 
i.e. the invariant mass of the dilepton pair. 

The generic letters 8 and indicate all the polar and 
azimuthal angles used in this work. The text and the figure 
captions specify each time which angles I am speaking 
about (suffixes resulted in difficult to read figure captions). 

Results will be presented with respect to angles mea- 
sured in two reference frames: 



a) "hadron collision frame" : the center of mass frame 
of the colliding proton and antiproton, with z along the 
proton direction. The other two axes are fixed once and for 
all, so this frame is the same for all events. The fragment 
momenta are always referred to this frame. The lepton 
momenta are referred to this frame when I look for corre- 
lations between them and the hadron momenta. 

b) "co-oriented dilepton cm. frame" : this is an event- 
dependent frame with all axes parallel to the former one, 
but boosted in such a way to be a center of mass frame 
for the lepton- antilepton pair. This frame is used to cal- 
culate the azimuthal cos (20)— asymmetry of the leptons, 
in a way that is as much similar to the tradition as pos- 
sible. In the present case, the main reason to avoid the 
Collins-Soper 23 frame and of other similar frames with 
event-dependent orientation of the axes is the difficulty in 
relating lepton and hadron properties. 

Another frame appears in the intermediate stages of 
the calculation: 

c) "spin frame" : an event-dependent frame, whose ori- 
gin and z— axis coincide with those of the hadron collision 
frame, and whose x and y axes are random-rotated with 
flat distribution around the z— axis. The quark and anti- 
quark spins are projected along the y— axis of this frame, 
but no initial or final observable quantity refers to it. 

2 General ideas 

I have built a "Drell-Yan-bascd" MonteCarlo model: as 
much as possible of the knowledge of the inclusive Drell- 
Yan process at partonic level has been used here. This 
includes the leading twist factorization properties [24 ,25, 
26,27 , the £— distributions and kr— distributions as they 
have been measured in this process and in the required 
regime [1611 17) . the phenomenology and theory of the in- 
clusive lepton cos (20) asymmetry [HIMISlIIIlinillSEllIini 
l2"gll29U30] . This should not hide the fact that I want to re- 
produce something that is not exactly a Drell-Yan process. 
In the individual channels composing the inclusive Drell- 
Yan process, higher twist effects may be much stronger 
than in the overall sum, and factorization properties have 
never been discussed. 

In the hadron-quark-spectator vertex I assume a ba- 
sic diquark form for the spectator, according with the 
widespread quark-diquark model used in [3T]. I will only 
consider one quark flavor (according to [20 90 % of the 
proton-antiproton Drell-Yan events derive from uu anni- 
hilations, so this is a good approximation). 

For the hadronization process, I will consider two al- 
ternative possibilities: (i) I will assume that spectator mo- 
menta may be identified with final hadron momenta, (ii) 
I will assume that the final qq pair, needed to allow the 
spectators to form physical hadrons, is created with a rea- 
sonably distributed relative momentum. 

My main goal is to reproduce properties of the final 
proton-antiproton pair, that have the same physical origin 
of the dilepton azimuthal cos(20)— asymmetry. 

In the pioneering works by [HUH] the cos (20)— asymmetry 
is one of the terms appearing in the most general form 
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for the dilepton pair distribution in unpolarized Drell- 
Yan. Within the scheme given in [28], this asymmetry is 
rewritten according with the quark-parton model, where 
it is proportional to the product of two terms (one for the 
quark, one for the antiquark) of the form h(x, kr)sT A kr, 
where h(x, kr) is the Boer-Mulders function. A commonly 
agreed definition for it may be found in [32], and some 
models exist |501|551|531|551|gg] . 

The qq — > 7* coupling selects qq pairs with opposite he- 
licitics. Combined with the Boer-Mulders spin-momentum 
correlation, this selects qq pairs with nontrivial space dis- 
tribution properties. The virtual photon is able to transfer 
some of these properties to the dilepton relative momen- 
tum, resulting in the known cos (2(/>)— asymmetry. 

To understand what this could mean for the spec- 
tator distribution, I trace this path back to its origin: 
the cos (2<^)— asymmetry derives from a property of the 
splitting vertex according to which a quark has nonzero 

< St A kr >• The spectators must balance this with 

< st f\pT,s P ect > — — < st A kr >, since in the splitting 
vertex PT, spec t + k T = PT, P arent = 0. This should imply 
azimuthal asymmetries of the final proton and antipro- 
ton, for the subset of events where these are observed in 
coincidence with a dilepton pair. 

The spin-momentum correlation picture of [28J poses 
some fundamental problems that have been solved at the- 
oretical level [30]. but there is some ambiguity about how 
to implement this theory in a parton MontcCarlo model. 
To clarify this, some general premise is useful: 

(i) In a theoretical model for a quantum process the 
expectation values of some variables are not obliged to 
have a physical meaning in the intermediate steps. In a 
MonteCarlo they should better have, since all the steps 
are classical, (ii) Loop integrals often present relevant fi- 
nal cancellations. In a MonteCarlo a single event is only 
one of the infinite configurations entering a loop. Some- 
times cancellations take place anyway over a large set of 
events, sometimes they do not. (iii) In a theoretical calcu- 
lation much is gauge artifact. In a MonteCarlo, only physi- 
cal particles and interactions can appear, (iv) Interference 
processes do not admit a straightforward implementation, 
(v) A "black box" distribution that turns around prob- 
lems always exists, but one would like to avoid it, as much 
as possible. 

The product of two spin-momentum correlations, each 
one implying < s y k x > 7^ for the quark and the an- 
tiquark respectively, does not violate general invariance 
laws. However, each of the two is individually unphysi- 
cal. In [30] the ultimate reason for the Boer-Mulders spin- 
momentum correlation is searched in initial state inter- 
actions between active and spectator partons originating 
from different parent hadrons. A nontrivial reduction work 
shows that it is possible to rewrite these effects as a distor- 
tion of the quark and antiquark initial-state momentum 
distributions, in a factorized format. 

This suggests two possible schemes for a MonteCarlo, 
that I name "strict factorization scheme" and "rescatter- 
ing scheme" . In the former one implements the conclusion 
of [30] 1 in the latter the starting point of the same work. In 



other words, in one case one assumes that the stationary 
momentum distribution of a quark in a hadron presents 
unphysical spin-momentum correlations, in the other case 
that this is the effect of interactions with partons coming 
from another hadron. 

At leading twist and given some constraints on the 
rescattering interaction, the two schemes lead to the same 
results (for the observable variables discussed in this work), 
for the following two reasons: 

(i) Since the works on leading-twist factorization in 
Drell-Yan 24,25,26,27 it is known that gluon exchanges in 
initial state interactions are in large part gauge artifacts, 
or cancel in loops, or may be reabsorbed in the distribution 
functions, or belong to a class of phenomena that have no 
effect in the limit of high collision energy. Translated into 
an implementation of rescattering, this means that only 
exchanges of transverse momentum at a fixed scale are 
admitted (a few GeV/c at most). This also implies kr -C 
k for all the partons. 

(ii) At leading twist, the qq — > 7* — > cross sec- 
tion (with assigned transverse polarizations for q and q) 
is symmetric for quark-antiquark exchange if the quark 
momenta respect kT <C k (see eq. [0J and the discussion 
in the Appendix). As a consequence, an event where an 
amount of transverse momentum is exchanged between a 
spectator and its companion quark, and an event where 
this momentum is exchanged between the same spectator 
and the oppositely coming antiquark, lead to the same 
observable final state variables. 

I will work in the rescattering scheme. It allows for 
organizing a chain of probabilistic steps, none of which is 
a priori unphysical. 

Another problem is posed by the fact that rescattering 
is the effect of an interference process. The probability of 
the considered process may be written in the form 

Prob ~ \(S*P*C + S*A*C)\ 2 . (1) 

Here 5" is the amplitude for the splitting process of the 
two hadrons and may be further factorized as S(1)S(2), P 
describes the free propagation of all the partons, A is first 
order rescattering, C is the hard qq — > l + l~ process, and 
"*" indicates sum/convolution with respect to the degrees 
of freedom of the partons in the intermediate stages. 

My assumption is that for any given spin configuration 
of the qq pair in the intermediate state I may rewrite 

\(S*P*C + S*A*C)\ 2 -> \S\ 2 * \P + A\ 2 * \C\ 2 

= |S(l)| 2 |S(2)| 2 *|^ + ^| 2 *|Cf(2) 

and substitute \P + A\ 2 with a physically allowed scatter- 
ing process between physical particles. 

An interaction that is suitable for this aim is the "scalar 
+ spin-orbit" one, well known in hadron scattering (see 
e.g. [37]). Its effect in the scattering of a spin-1/2 projec- 
tile by an unpolarized target is to produce a final mo- 
mentum distribution of the form a(fc^) + b(kT)sT A fcr- It 
does not violate any invariance law, derives from an inter- 
ference process between scattering and no-scattering am- 
plitudes and modifies the quark momentum distribution 
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in the same way as an "intrinsic" Boer-Mulders function 
would do. 

What is done here is the probabilistic version of the 
quantum method used in [JO], where one rewrites at par- 
ton level the hadron-hadron spin-orbit interactions. The 
method has some overlap with the 1-gluon exchange used 
in perturbative models 30,33,34,35,36 . 

Alternatively, spin-orbit rescattering may be just con- 
sidered an effective way to distort the quark momentum 
distribution, so to introduce a Boer-Mulders term oc Sx A 
kx with pre-assigned shape, without worrying about its 
physical explanation in more fundamental terms. 



3 Constraints from the features of inclusive 
Drell-Yan 



The simulation code is parameterized in such a way 
that the final outcome reproduces the features of eq[5] and 
of the related phenomenology for the dilepton distribu- 
tions in n~W fixed target experiment s [ 1 Oil 1 1 II 1 2 II 1 3 II 1 4) . 
We have not data on the cos(2(j>) asymmetry in pp Drell- 
Yan, so I have to rely on tt~ W data. In both cases we have 
dominance of valence- valence uu annihilations. 



4 MonteCarlo structure 

The event are sorted at s = 100 GeV 2 . Quarks and lep- 
tons are considered light-like. With regard to spectators, 
we only need their 3-momentum, so it is not relevant to 
specify their mass. All the events are subject to the cutoff 
r > 4 2 /100, where tekk M 2 /s. 



A detailed discussion of the following equations may be 
found in a very extensive form in e.g. [38] . and in more con- 
cise and specific form in |14II29| . or in many other works. 

Neglecting quark spin effects and assuming one rel- 
evant quark flavor, the parton model PP —> + X 
cross section is 



A 



u(x, k T )u(x, k T ) B. ta ,V lv 



(3) 



Here A is a normalization constant, the factor 1/xxs « 
1/q 2 combines flux and phase space factors, u[x, kx) is 
the u — quark unpolarized distribution in an unpolarized 
proton. L^ v is the lepton tensor, is the qq tensor, 
defined so that H^L^ V / q A is adimensional. 

The quark and antiquark variables entering this equa- 
tion are x, x, kx and kx- From the measurement of the 
lepton momenta p and p we may reconstruct x, x and 
qx = Pt + Pt = kx + kx- The individual quark trans- 
verse momenta are not observable, but some features of 
the distribution of kx — kx may be inferred from the 
distribution of px — pr- 

To introduce quark spin effects in the generator I rewrite 
the above cross section in the form 



u(x, kx)u(x, k T ) H^Sy, s y )L^ v 



.(4) 

where a random axis y is chosen, and the polarizations 
s y , s y are random sorted with distribution /(%). 

Eq|4] gives the same results as eq|3l as far as the mo- 
mentum and the spin distributions u(x, kx) and f(s y ) are 
uncorrelated. The spin-orbit rescattering, or equivalently 
the Boer-Mulders terms, introduces a correlation. In this 
case the phenomenological distribution is 



da 



l + cos 2 {9) 



xxs 



(u{x, kx)u(x, kx) 



h(x,k T )h(x,k T )sin 2 (0)cos(2(t))) d[cos{9)]dcj>. (5) 
(see sections 4.4 and 5.1 for the definition of and (f>). 



4.1 General flowchart 

The chain of steps in the simulation is: 

1) The proton with initial momentum Pqz w y/s/2z 
splits into a light-like quark and a diquark. The quark 
has 3-momentum [k x , k y , xPq] that is sorted randomly: x 
is sorted according to a distribution f(x)/x; k x and k y 
refer to the random-oriented spin frame, and are sorted 
with gaussian distribution with center at k x — k y = 0. 
Details on the parameters are given later. The factor 1/x 
in f(x)/x is needed to take into account a factor 1/xx in 
the cross section (eqJ3J). 

2) The 3-momentum of the spectator diquark is set to 

[~kx, ~ky, (1 - X)P }. 

3) A value s y = ± 1 (meaning y— spin = ±h/2) is 
sorted for the quark, with probability 50% vs 50 %. At this 
stage, the sorting of the spin and of the other variables are 
independent. The diquark is assumed a scalar one, so it 
has no spin. As far as final state rescattering is neglected, 
the difference between scalar and vector diquarks is not 
relevant. 

4) The same operations (1, 2, 3) are performed for the 
antiproton, whose initial momentum is — Pqz. 

5) Spin-orbit shift of the quark momentum: k x — > k x + 
s y Ak x . The shift Ak x is gaussian-distributed, with posi- 
tive average value of < Ak x >. After this step the quark 
has < k x > — s y < Ak x >. 

6) The anti-diquark (i.e. the spectator coming from 
the antiproton) is subject to the opposite shift: its ini- 
tial transverse momentum kd, x = —k x is shifted to —k x — 
s y Ak x . The result of this and of the previous step is that 
the total momentum of the quark and of the anti-diquark 
is conserved. 

7) Steps (5,6) are symmetrically performed on the q- 
diquark pair. 

8) The quark and the antiquark variables including 
spins are used to build the hadron tensor. 

9) In the lepton co-oriented cm. frame (see section 
1.1), for q 2 fixed by the quark and antiquark momenta, 
the angles specifying the direction of the lepton momenta 
are sorted isotropically. After transforming them to the 
spin frame, the lepton tensor is build. 
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10) The lepton and the hadron tensor are contracted, 
and the scalar W(s y ,s y ) = H^ v ' (s yi s y )L au / 'q 4 is used in 
an accept/reject procedure, where the event "qq trans- 
formed into is accepted or rejected with probability 
oc W(s y ,s y ) (W does not depend on q 2 , but only on the 
relative orientation of the momenta). If the event is re- 
jected, one restarts from step (1) and all the variables 
have to be sorted again. 

11) In the case of an accepted event, the observable 
momenta (still in the spin frame) are transformed to the 
hadron collision frame or to the lepton co-oriented cm. 
frame for the data analysis. 



4.2 fey-distributions and s x k y distortion 

Here I discuss the sorting of the quark variables. For the 
antiquark, the procedure is the same. 

The initial components of kx are sorted with gaussian 
distribution, with < k x > = < k y > = 0, and y<fe|T> 
= 0.4 GeV/c. 

The shift Ak x is sorted randomly with quasi-gaussian 
distribution with < Ak x > = 0.5 GeV/c and fluctuation 
0.27 GeV/c. "Quasi-gaussian" means that we have the 
boundary conditions —1 GeV/c < Ak x < 2 GeV/c, and 
that the sorting probability tends regularly to zero at these 
limits. The value 0.27 GeV/c has the only motivation of 
reproducing the phenomenological data of [T0l[TTl [T2, 13, 

eng. 

The spin-orbit effect for quarks from a proton and an- 
tiquarks from an antiproton is here assumed to have the 



same sign of < (st A kx) ■ P 



parent 



>. This means that 



< (st A kx) ■ z > is opposite in the two cases, since the 
two parent hadrons are opposite-directed. 



4.3 a; distribution 

The a;— distribution has the basic form f(x) = x a (l — x) b . 
Fixed-target Drell-Yan experiments at s ~ 100-600 GeV 2 
used the same parameterization and gave a ss 0.5, 6 ~ 2 
(see e.g. the appendix of [14]). I have adopted these two 
values. 

It must be noted that the lower limit sxx > (4 GeV) 2 
weakens the experimental constraints on a, and decreases 
the relevance of a precise value for it. 



4.4 Polarized quark lepton tensor contraction 

The calculation of this tensor-tensor contraction is re- 
ported in the Appendix. For steps (8,9,10) of the previous 
flowchart I need the probability distribution 

W{Sy,Sy) = H^{ Syi S y )L^/q i (6) 



where L av is the lepton tensor and H^(s yi s y ) is the ten- 
sor associated to a quark and an antiquark with assigned 
transverse polarizations. 

Extracting a factor q from each vector: 



I get 



and 



IV, 



unpol = g [(Nn)(Nn) + (Nn){Nn)}. 



(8) 



= [1 + (SySy)} ■ W, 



/ unpol 



s y Sy(l + n y ny) (9) 



with s y = ±1. To get to this result, terms with magni- 
tude kr/q have been systematically neglected. So eq|H]is 
a leading twist result. 

The first term in eqE] is about 3 times more relevant 
than the second one in determining the overall cross sec- 
tion. In the center of mass of the qq — > annihilation, 
where the virtual photon is at rest, and the 3-vectors n, 
AT, etc, are unitary, the unpolarized term W unpo i favors 
configurations where the quark and the lepton directions 
are aligned, according to the known "1 + cos 2 (0)" law. 

In fig.l I show the distribution f(\cos(9)\) together 
with a fitting curve oc [1 + cos 2 (9)} where 9 is the lepton 
polar angle in the lepton co-oriented cm. frame. In the 
Collins-Soper frame we should get a precise [l+cos 2 (9)]— law, 
and for large s the z— axes of two frames coincide. Fig.l 
shows that we are very close to this condition. 




0.8 0.9 1 

cos(6)| 

Fig. 1. Event distribution w.r.t. the polar angle of the positive 
lepton in the lepton co-oriented cm. frame. The continuous 
curve is a fit of the form iV = No[l + cos(6) 2 ]. For each point, 
a very small error bar is present although not visible, estimated 
by the standard fluctuation. 



k^ = (q/2)N fl 1 



(«7/2K, 



5 Results 

I evaluate the measurable asymmetries according to the 
definition 

_ N(+)-N(-) 
A± = N(+) + N(-) (10) 

where N(±) is the simulated number of events correspond- 
(7) ing to a positive or negative value of some observable. 
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In most of the following cases the observable is a factor 
cos(2cj>), where (f> is a lepton or hadron azimuthal angle 
ranging from to 2-7T. I also present figures with event 
distributions vs 6, and hand-made fits of the form 



events 



f{N,P,<f>) = N l + /3cos(2<j>) 



(11) 



For a distribution exactly coinciding with eq llll the asym- 
metry calculated according with eqUU] is 



IT 



(12) 



A typical value that I find is (3 rj 0.2 corresponding to A± 
« 0.13. 

All the following simulations have been performed once 
more after removing any spin-momentum correlation. Be- 
cause of statistical fluctuation, or of uncontrollable nu- 
meric systematic effects, these "vacuum" simulations show 
asymmetries of magnitude 2 %. So we know that few- 
percent magnitudes cannot be taken seriously given the 
event numbers used in the following. 



5.1 Behavior of the lepton cos{24>)— asymmetry 

In this subsection <j) is the angle between the projections 
on the transverse xy plane of the vector q = p + p mea- 
sured in the hadron cm. frame, and the vector p — p mea- 
sured in the lepton co-oriented cm. frame (where it simply 
coincides with 2p) . This is equivalent to the ordinary def- 
inition of the angle appearing in the cos {2(f>) — asymmetry, 
although the traditional procedure is to rotate the lepton 
cm. frame in such a way that <f> is just the angle be- 
tween the lepton transverse momentum and the x— axis. 
Some slightly different choices are normally adopted for 
the z— axis, whose implications have been discussed in |39j . 
In [Tl] the data have been organized according to different 
frames, and one may appreciate the differences. 

In fig. 2 I show the (^—distribution for the events with 
qx between 1.5 GeV/c and 3 GeV/c. The fitting curve is 
1 + 0.23 * cos(2(f) (apart for an overall normalization). 

The applied rescattering model does not introduce di- 
rectly any x— dependence or qx— dependence: all the quarks 
are (statistically) subject to the same spin-dependent mo- 
mentum shift. However, indirectly a strong qx— dependence 
and some x— dependence is produced in the final results. 

In Table 1 I report the event numbers and the asym- 
metry of a set of simulations corresponding to a series of 
consecutive bins in qx- The event numbers in this table are 
proportional to the corresponding cross sections. For these 
event sets a further cutoff has been applied on the polar 
angle: only events with \9\ < 45° have been counted in this 
statistics. Apart for enhancing the cos (2cf>) — asymmetry, 
the second choice is useful to avoid the problems described 
in [41], and similar cutoffs have been adopted largely in 
the previously quoted experiments because of small-angle 
and large-angle acceptance problems. 
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<I> (71 units) 



Fig. 2. (^—distribution of the simulated dilepton events, with 
the transverse momentum of the virtual photon satisfying 1.5 
GeV/c < qr < 3 GeV/c. For the definition of (j) see the text. 
The fit is oc 1 + 0.23 * cos(2(j>). 



Range 


Event 


Event 


Asymmetry 


(GeV/c) 


difference 


sum 




0-0.5 


163 


32564 


0.005 


0.5-1 


1093 


56314 


0.02 


1-1.5 


2035 


33826 


0.06 


1.5-2 


2061 


12677 


0.16 


2-2.5 


872 


3435 


0.25 


2.5-3 


208 


616 


0.34 



Table 1: Distribution w.r.t. qx of dilepton events 
and of their cos (20)— asymmetry, with the general 
cutoff \tan(9)\ < 1. For the definition of the 0— angle, 
see the text. 



The chosen parameters of the distorting term produce 
a cos {2(f)— asymmetry with a reasonable magnitude and 
qx— dependence (see the data in pT0irmfT2lfT4] . and the fit 
in [25]), at least up to 2.5 GeV/c. 

The experiments [OTfmn^lfTl"] show a decreasing trend 
at large x, and give no information for x below the valence 
region. The proton-proton measurement by E866[T5] had 
larger beam energy and smaller average x and x, and in- 
volved one sea active quark at least. It showed near-zero 
asymmetries. 

My simulations show a decrease of the asymmetries for 
large and small x, x. It is however difficult to disentangle 
the role of the two. To properly analyze the small— a; case 
I should remove the limit q > 4 GeV/c 2 on the dilepton 
mass. This would take me into the region qx/q ~ 1, where 
I cannot rely on some approximations used in building the 
simulation code. Also, near the edges of the phase space 
the physics may be seriously affected by constraints that 
I have not applied here, and by higher twist terms. So I 
would not take seriously simulations in the regions where 
x and/or x, are close to zero or one. In the other regions, 
the asymmetries are roughly x— independent. 
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5.2 Proton-antiproton cos {2<f>) asymmetry 

In Table 2, I report final hadron asymmetries and com- 
pare them with the corresponding lepton asymmetries in 
two relevant kinematic situations. For the hadron cos(2(f>) 
asymmetry, <f> is the angle between the difference and the 
sum of the transverse momenta of the proton and of the 
antiproton, in the hadron collision cm. frame (remark: 
the sum is — qx)- For the lepton cos (20)— asymmetry, the 
angle is the same as in the previous section. Apart for the 
transverse momentum cuts, and for the cutoff 4 GeV/c 2 
on the dilepton mass, no further cuts have been applied 
to these simulations. 



Pair 


Range 


Event 


Event 


Asymmetry 




(GeV/c) 


diff. 


sum 




leptons 


1.5 - 3 


17921 


133498 


0.1342 


hadrons 




17847 


133498 


0.1336 


leptons 


- 1.5 


16763 


974911 


0.017 


hadrons 




48550 


974911 


0.050 



Table 2: Compared cos(2(j>)— asymmetries for lep- 
tons and hadrons 



For the first case of Table 2, the 4>— distribution of 
the hadron pairs is reported in fig. 3. The fitting curve is 
1 + 0.22 * cos(2(f>) (apart for an overall normalization). 
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Fig. 3. Distribution of the proton-antiproton pairs w.r.t. the 
angle <f> between the sum and difference of the hadron trans- 
verse momenta in the hadron collision frame. These events are 
the same as in Table 2, upper part, with cutoff 1.5-3 GeV/c for 
the total transverse momentum and 4 GeV/c 2 for the dilepton 
mass. The fit is oc 1 + 0.22 * cos(20). 



5.3 hadron-lepton correlations 

At leading twist the hadron and the lepton vectors can- 
not communicate through quantities that are odd in the 
charge exchange on one of the two sides (hadron or lep- 
ton). I have checked the correlations (p x py—PyPx)iept(PxPy— 



P y Px)hadr, and (p-p)i ep t ■ (p-p)hadr (transverse compo- 
nents), for which I have found nothing meaningful. 

The 0— angle between the relative transverse momen- 
tum of the leptons and of the hadrons could present (f>— even 
systematic behaviors, since it is the difference between 
4>i ep t and Aspect when both are measured in the same 
frame. However, the 4>— distribution in fig. 4 does not show 
reliable systematic deviation from homogeneity. I remark 
that in fig. 4 the vertical event scale starts from 5000. If 
it started from zero as in the other figures, it would ap- 
pear almost completely flat. The fitting curve is 1 + 0.03 * 
cos (2(f>) (apart for an overall normalization). This corre- 
sponds to a possible asymmetry about 2 %. This mag- 
nitude is within the simulation uncertainties and is not 
meaningful. 

events 1 



7000 



6000 7 
5500 
5000 - 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 

3> (71 units) 

Fig. 4. Event distribution w.r.t. the angle 4> between the rel- 
ative transverse momentum of the diproton pair and the rela- 
tive transverse momentum of the dilepton pair, both referred 
to the hadron collision frame. Same cuts as in fig. 3. The fit is 
oc l + O.O3*cos(20). 



5.4 Single hadron - leptons correlations 

A strong correlation is the one between one of the hadron 
fragments (e.g. the proton) and the lepton momentum (p— 
p)iept (transverse components). 

From a theoretical point of view this adds nothing new, 
because the momentum of any of the spectator partons 
may be written as — q/2 ± (p — p) spect /2, so any observ- 
able that is linear with respect to the momentum of a 
spectator is not independent from those ones that I have 
considered previously. From the experimenter's point of 
view, however, to measure a proton-dilepton asymmetry 
in the laboratory frame is easier than detecting both a 
proton and an antiproton (the leptons must be detected 
in any case). According with the analysis in [50], in the 
specific case of PANDA we have a large number of events 
where it is possible to identify a \x~pp— event from the 
detection of /z + , fir and p only. 



8 



A. Bianconi: Azimuthal asymmetries in exclusive four-particle "Drell-Yan" events. 



In fig. 5 I report the distribution of the 0— angle be- 
tween the proton momentum and the difference between 
the lepton momenta (all in the hadronic center of mass 
frame). The only cutoff is qr > 1.5 GeV/c. The corre- 
sponding cos(2(j))— asymmetry is 0.14 to be compared to 
the value of the lepton asymmetry that is 0.13. The fit is 
oc 1 + O.22*cos(20). 
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Fig. 5. Event distribution w.r.t. the angle <fr between the pro- 
ton transverse momentum and the dilepton relative momentum 
in the hadron collision frame. Same cuts as in figs. 3 and 4. As 
remarked in the text, the evident asymmetry is just another 
way to measure the ordinary dilepton cos(2<f>)— asymmetry. 
The fit is oc 1 + 0.22 * cos(2cf>). 



6 Higher twist and hadronization effects 

A very simple setup for a realistic hadronization is the one 
used by Pythia in the small phase-space limit of the string 
fragmentation model: the creation of a single qq— pair with 
a random relative transverse momentum (see the detailed 
description in section 2.4.1 of |4"2"] . and eq.25 in [13]). I 
have simulated the effect of this on the distribution of the 
final hadron pair. 



Gaussian width of 
the relative Pt due 
to hadronization 


Asymmetry 





0.134 


0.35 GeV/c 


0.121 


0.7 GeV/c 


0.091 



Table 3: Changes of the hadron cos (2(f)— asymmetry 
associated with a relative gaussian-distributed trans- 
verse momentum, due to the qq pair that forms the 
final hadrons together with the diquarks coming 
from the initial state. 0.134 is the reference value 
from Table 2. 

In other words, my previous simulation relative to table 
2 (upper line) has been repeated with the addition of a 



gaussian broadening of the relative transverse momentum 
of the final hadron pair. I have tested two cases, where the 
gaussian broadening has width 0.35 and 0.7 GeV/c. The 
results for the hadron cos(20)-asymmetry are reported in 
table 3. 

Since 0.2-0.4 GeV/c is a reasonable range for the width, 
we see that the hadronization process, at least in this sim- 
ple form, does not change the hadron azimuthal asymme- 
try. 

Some authors (see 44,45,46 ) have proposed modifi- 
cations of this basic form of the string model, so to in- 
clude spin-dependent mechanisms in a 2-hadron produc- 
tion stage. Probably, this class of processes would more 
seriously affect the hadron asymmetries calculated by me, 
but they cannot touch the lepton asymmetries. So, the 
approximate equality of lepton and hadron asymmetries, 
predicted in this work, is a signature of the absence of 
final state spin-dependent effects. 

The process pp — > +pp has chances to be affected 

by higher-twist effects, since it is an exclusive process. 
In a scheme where at least approximate factorization is 
present and higher twist terms are a correction and not 
the main part of the cross section, these effects may be 
separated into initial-state interactions (that link active 
and spectator partons) and final-state interactions (not 
implying direct or indirect momentum exchange between 
spectator and active partons). 

Higher-twist effects in the initial state should preserve 
the approximate equality of hadron and lepton asymme- 
tries (although the shape of both could differ from what I 
suggest here). Indeed, the momentum-conservation mech- 
anism that is behind the dilepton-dihadron correlation is 
independent of the details of the initial-state processes. 

Summarizing: lepton cos (2</>)-asymmetries that are in 
disagreement with my predictions (fig. 2, table 1) are a 
signature of higher twist effects in general, while the in- 
equality of hadron and lepton asymmetries is a signature 
of specific spin- dependent final state effects. 

7 Conclusions 

This work has considered exclusive pp — > fi + ji~ +pp events 
at s = 100 GeV 2 . The main focus is (i) the possibility of 
a lepton cos(2(f>) — asymmetry, of the same nature of the 
one observed in inclusive \x + \x~ + X events at the same 
energy, (ii) the possibility of associated asymmetries of the 
recoiling hadrons. 

The model simulations presented here show that cos (20) 
asymmetries of the leptons in their cm. frame should be 
accompanied by cos(2<f) asymmetries of the hadrons in 
the laboratory frame. In the latter case, </> is the angle be- 
tween the difference and the sum of the hadron transverse 
momenta. The size of these asymmetries is the same. 

If the measured asymmetries differ from what I have 
suggested here, this is a signature of higher twist effects 
(since my predictions for the lepton distributions may be 
read as a re-parameterization of the known distributions 
measured in inclusive Drell-Yan). But, if the lepton and 
hadron asymmetries differ strongly between them, this is 
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a more specific signature of final state spin-dependent ef- 
fects. 



A Appendix: Polarized quark— lepton 
contraction H llv (sy,Sy)L^. 



We need the contraction 



q 4 W 



H^L„ 



(13) 



of the quark-level hadronic and lepton tensors in the pro- 
cess qq —> , where all four particles are ultrarelativis- 
tic. I use the shortened notation for the traces 



T[..] 



\Tr[..]. 



(14) 



I use the definitions of the Berestevskij-Lifsits-Pitaevskij 
book[47] (better known as the 4th book of the Landau- 
Lifsitz Course in Theoretical Physics) . As the only excep- 
tion to this, I indicate k^j^ with the widespread notation 
# instead of using k as was done in that book. 



A.l The case of unpolarized quarks 

If the quarks are unpolarized, we simply get 

q 4 W unpol = T[# 7 ^ 7 H T^^] = 

= 2 [{kp){kp) + (kp)(kp)}- (15) 

I extract a factor q from each vector: 

k» = (q/2)N", p» = (g/2)n", ... (16) 

getting 

W unpol = I [(Nn)(Nn) + (Nn)(Nn)}. (17) 
o 

In a center of mass frame of the partonic process, the 
3- vectors TV, n etc are unitary vectors. In the hadron 
collision frame this property is not true, but of course 
the full tensor contraction is frame-independent so W unpo i 
is a scale-independent object. The same is valid for the 
leading-twist terms in the spin-dependent W(s y ,s y ) that 
will be calculated below. 



axis that accompanies each individual quark in its move- 
ments, so I have assigned a common quantization axis (the 
y— axis) for all spins. Since here quarks present nonzero kx 
components, this means 0(kx/q) helicities. 

The relation between s M and the polarization in a rest 
frame <x is, for massive particles, 

so = cr L \k\/m, s L = a L E/m, s T = cr T (19) 

and evidently it creates problems for U.R. particles, un- 
less the longitudinal component is strictly zero. However, 
when the previous expressions are used to write the den- 
sity matrix eq |18l in terms of the rest frame polarizations, 
E/m— terms cancel and we are left free from mass singu- 
larities: 



#(l- 7 5 ^) = #[l-j 5 (±a L +* T -y T )} 



(20) 



(± differentiates particles and antiparticles) . This says 
that a way exists to make calculations without meeting 
mass singularities. 

In the following it will be comfortable to use invariant 
trace formalism. So, it is useful to understand which ap- 
parently large terms may appear and how to treat them. 

In the product of two (1 — 7 5 />) density matrices with 
not exactly transverse spins, special care is needed with 
the terms containing the product of one or two longitudi- 
nal components. The former give zero trace for odd par- 
ity of the number of gamma matrixes. The latter modify 
W unpo i to 



W h ~ h = (l-hh)W unpol . 



(21) 



where h and h are the average helicities due to the longi- 
tudinal projections of the quasi-transverse spin. The mag- 
nitude of these components is 



h ~ ky/q, h ~ ky/q, 



(22) 



that leads to W h ^ = W unpo i + higher twist terms, sys- 
tematically neglected here. 

The full trace for polarized quarks is 



H: 



[IV 



7W(l- 7 5 ^) P7 1/ (1-7 S 



(23) 
(24) 



This gives a sum of 15 terms. Two of them contain the 
products (sk) = and (sk) = 0. We are left with 



A. 2 polarized quarks (leading twist) 

In the case of a polarized (anti)quark 
# -> #(l-7 5 ^) 



(18) 



+ (fcs){/P,<f} + {ks){k",s v } - (ks)(ks)g flv (25) 



where s M is the polarization 4-vector for the quark, re- where {a^frj,} = a u b v + a v b p 



specting the exact 4-dimensional constraint (sk) = 0. If 
the spin is strictly transverse to the particle motion, then 
also s ■ k = 0. 

In the implementation of a multi-particle MonteCarlo, 
it would not be practical to define a spin quantization 



After contracting the hadron tensor with the lepton 
tensor, the final result consists of many terms. I drop the 
hh— terms on the ground of eq. 1221 The transverse spin 
dependent terms contain the product q s y s y , and have 
forms like e.g. q A s y s y N y n y (Nn). I neglect systematically 
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all the products that contain one at least among N x , N y , 
N x , N y . These terms do not reach the magnitude q A , since 
N x ~ k x /q etc. 

I retain all the terms containing (ns){ns), (ns)(ns), 
(ss) (transverse components of n and n with magnitude 
~ 1 arc statistically frequent). 

The final result is 

W Sv a- v = [l + (SySy)]-W un pol - 7;SySy(l + n y n y ). (26) 
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